Nambu-Goldstone Mode in a Rotating Dilute Bose-Einstein Condensate 
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The Nambu-Goldstone mode, associated with vortex nucleation in a harmonically confined, two- 
dimensional dilute Bose-Einstein condensate, is identified with the lowest-lying envelope of octupole- 
mode branches, which are separated from each other by the admixture of quadrupolar excitations. 
As the vortex approaches the center of the condensate and the system's axisymmetry is restored, 
the Nambu-Goldstone mode becomes massive due to its coupling to higher rotational bands. 
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One of the unique features of a gaseous Bose-Einstein 
condensate (BEC) is the fine tunability of interatomic 
interactions due to the Feshbach resonance [1, 2]. This 
degree of freedom can be utilized to explore some unique 
features of a rotating BEC by making the strength 
of interaction close to zero; then as the angular mo- 
mentum (AM), L, of the system increases, the low- 
lying states of the BEC in a harmonic potential become 
quasi-degenerate [3, 4] and hence highly susceptible to 
symmetry-breaking perturbations. Such high suscepti- 
bility is considered to be the origin of vortex nucleation. 
Both experiments [5] and mean-field theories [6, 7] have 
demonstrated that as L increases, vortices enter the sys- 
tem from the outskirts of the BEC by spontaneously 
breaking the system's axisymmetry. Associated with this 
symmetry breaking must be a Nambu-Goldstone mode 
(NGM) which, however, has been elusive. This paper re- 
veals the NGM using many-body theory and shows that 
this mode becomes massive as the vortex approaches the 
center of BEC at which point the axisymmetry of the 
system is restored. 

We consider a system of N identical bosons, each 
with mass M, that undergo contact interactions 
and are confined in a two-dimensional harmonic 
potential with frequency u>. Throughout this pa- 
per, we measure the length, energy, and AM in 
units of (h/Mu) 1 ' 2 , hu>, and h, respectively. The 
single-particle and interaction Hamiltonians are then 
given by H = £f =1 {-2d 2 /d Zj dz* + \ Zj \ 2 /2) and 
V = 2TtgJ2j=£kd( z j ~ z k), where Zj = Xj + iyj is the 
complex coordinate of the j-th particle and g gives the 
ratio of the mean-field interaction energy per particle 
to hu. When g -C 1, it is legitimate to restrict the 
Hilbert space to that spanned by the basis functions 
(f) m (z) = (z m /\W)e-l ! l 1 / 2 (m = 0,1,2,-..), where 
m is the AM quantum number. In this "lowest- 
Landau-level" approximation, the field operator is 
expanded as \E , (z) = J2m=o bm 4> m {z), where b m is 
the annihilation operator of a boson with AM m. 
The second-quantized form of V then reads [8] V = 

9'^2im 1 ,---,m i Vrn 1 ,---,m 4 bl ni bl n2 b m3 b m4 , where V mx , — ,m,i = 



S mi + m2 ,m 3 + mi (mi + m2)\/(2 mi+m2 y/mi\m 2 \rn 3 \m 4 [). 
Since Hq = L + N is constant for given L and N, the 
dynamics of the system is determined by V alone. The 
lowest-energy state of this system at a given L is referred 
to as the yrast state, whose interaction energy is given 
for a repulsive case by gN(N — 1 — L/2) [8, 9]. The trace 
of the yrast state viewed as a function of L is called the 
yrast line. In the following we will measure the energy 
of the system from that of the yrast state. 

NGM in a rotating BEC-. The NGM should appear 
when the vortex is about to nucleate; that is, when 
the axisymmetry of the system is being broken. In this 
regime (L <C N), the excitation spectrum is divided into 
two groups having different characters [10]. One of them 
involves excitations whose energies are on the order of 
gN, and the other involves pairwise repulsive interac- 
tions between octupole modes with excitation energies 
given by 6 = (27s/34)SJS 3 (SJS 3 - 1) [10]. Because the 
energy scale (~ g) of the latter group is smaller, by a 
factor of 1/N, than that of the former one and hence 
vanishes in the thermodynamic limit, one might suspect 
that O is the NGM. However, this is not the case. We 
show that the NGM is the envelope of equally-spaced 
octupole-mode branches. The envelope is labeled G in 
Fig. 1, which is obtained by exact diagonalization of V 
for N — 256. This beautiful structure was not found 
previously because it emerges only for large N and for 
relatively large values of L/N [11]. In Fig. 1, octupole- 
mode branches are indicated by solid and dotted curves, 
which are equally spaced with AL = 2. We have con- 
firmed that this spacing is caused by the admixture of 
one quadrupole-mode excitation. 

For a mode to qualify as a NGM, it must meet three 
conditions. It must be massless, it must be associated 
with a broken symmetry, and it must play the role of 
restoring the broken symmetry. Because the G mode 
belongs to the second group discussed above, the ex- 
citation energy vanishes in the thermodynamic limit. 
Thus the first condition is met. Since the yrast states 
for L < N are energetically degenerate under a rota- 
tion with angular velocity fl = ^jjjf = lu(1 — 
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an axisymmetry breaking is expected to occur even for 
an infinitesimal, anisotropic symmetry-breaking pertur- 
bation. In fact, as the vortex enters the system, the 
axisymmetry is broken as seen in mean-field density 
profiles in Fig. 2 (a). Thus the second condition is 
met. The symmetry-restoring force is a correction to 
the mean-field result, so that many-body theory must 
be invoked to find out whether the third condition is 
met. We note that in exact diagonalization calculations 
the single-particle density-distribution function, p(r) = 
(^ t (r)^E'(r)), is isotropic [i.e., p(r) — p(\r\)] and is not 
suitable for studying the problem of axisymmetry break- 
ing. The symmetry breaking can be studied with the 
conditional distribution function (CDF) [12] defined as 
p(r;r ) = ^ (r)& (r )$ (r )$ {r)) , whe re r is the 
position of a test particle. In order to ensure that the 
back action of placing the test particle on the system is 
negligible, we take |ro| =3 which is well outside the con- 
densate. We set ro = (ro,0) without loss of generality. 
Expanding the CDF in terms of cos(m0), we have 

N 

pO; r o) = -t-j — rp(r;r )-p(|r|) 



= 2^ C m (r;\r \) -=—, (1) 

m=0 v 

where r and 9 are the polar coordinates of r. 

Figures 2(b) and (c) illustrate the CDF p(r; |r | = 3) 
for the lowest-lying excited states of 76 bosons with L = 6 
and L = 12, respectively. The CDF in Fig. 2(b) fea- 
tures three-fold symmetry, reflecting the fact that the 
AM of this state is carried mostly by octupole-mode ex- 
citations. We also note that one of the peaks is located 
near the boundary of the condensate where the vortex 
comes in. This implies that the octupole-mode excita- 
tions compensate for the density depletion caused by the 
entrance of the vortex and thus plays the role of restor- 
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FIG. 1: Nambu-Goldstone mode (labeled G) and its mass 
acquisition in a rotating BEC. The energy is measured from 
the yrast line (horizontal bottom line). The solid curve shows 
the branch that arises from octupole-mode excitations alone 
and is described by (27g/34)ri3(n3 — 1), where ri3 denotes the 
number of octupole-mode excitations. The solid and dotted 
curves are obtained from the least-squares fit of the data to 
quadratic polynomials. The dotted curves are displaced from 
the solid one by 2, 4, 6, 8 units of angular momentum. These 
shifts are caused by the admixture of quadrupole-mode exci- 
tations. The lowest-energy envelope of these branches consti- 
tutes the Nambu-Goldstone mode which acquires mass as L 
increases, as evident from the curve labeled G. The excited 
states that involve the center-of-mass motion are not shown. 



FIG. 2: (a) Mean-field density profiles of BEC for L/N = 
0.1,0.4, and 0.8. (b)-(c) Conditional distribution functions 
(CDFs) p(r; |ro| = 3) defined in Eq. (1) for the quasi- 
degenerate lowest-lying excited state with L = 6 (b) and 
L = 12 (c). (d) CDF for the yrast state at L = 12. (e)- 
(g) show expansion coefficients C m (r; \ro\ = 3) for m — 1, 2, 
and 3, respectively, for the lowest-lying excited state for each 
L satisfying 8 < L < N = 76. 
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ing the axisymmetry of the system. A further increase in 
L is caused by quadrupolc-mode excitations (see dotted 
curves in Fig. 1), so that the lowest-lying excited state 
gradually loses the character of three-fold symmetry, as 
illustrated in Fig. 2(c), and for larger L the CDF even- 
tually shows a dipole-like structure. On the other hand, 
Fig. 2(d) shows the CDF of the yrast state at L = 12, 
featuring a dipole-like structure. The dipole-like distri- 
bution supports the entrance of the first vortex into the 
condensate by breaking the axisymmetry of the system. 

Figures 2(e), (f), and (g) show the expansion coeffi- 
cients C m (r:\r \ = 3) in Eq. (1) for m = 1, 2, and 3, 
respectively, for the lowest-lying excited states, which 
are constituting the G branch in Fig. 1. We note that 
Ci(r; |r | = 3) is negative with large magnitude around 
r = 1 (i.e., around the periphery of the condensate), 
which suggests that this component is responsible for the 
invasion of the vortex. On the other hand, the signs of 
C m (r; | ro | = 3) for m = 2 and m — 3 are positive for 
L <C N, and their magnitudes are maximal around r = 1. 
This suggests that these two modes cooperate with each 
other to counteract the density depletion caused by the 
entrance of a vortex, in agreement with what is stated 
above. We note that the amplitudes of C2 and C3 dwin- 
dle rapidly and that Ci changes its sign when L increases 
beyond the quasi-degenerate region. This, too, supports 
our claim that they constitute the NGM. We have con- 
firmed that the signs of C4 and C5 are also positive, but 
that their magnitudes are negligible compared with that 
ofC 3 . 

Mass Acquisition of the NGM-. Finally, we investigate 
the mechanism by which the NGM becomes massive 
as L increases. To this end, we separate the coordi- 
nates into center-of-mass (COM) z c = (l/AOEjLi-Zj 

and relative ones Uj = Zj — z c , where Eili u i = 0- 
Accordingly, H is decomposed into the COM part 
H com = -(2/N)d 2 /dz c dz* + N\z c \ 2 /2 and the rela- 
tive part H rcl = (2/N)J2"-M l - N5 jk )d 2 /du du* k + 



E 



JV 



-/2. The total Hamiltonian is thus decomposed 



into the COM part H com and the relative part H rcl + V, 
where V = 2ngJ2j^k ^( u j ~ Uk )- Therefore, a many- 
body wave function is separated into their counterparts: 
y(z u z 2 ,---,z N ) = * com (2 c )* rel (ui,n 2 ,---,« JV -i). In 
what follows, we will ignore the COM part because it is 
decoupled from V. 

We introduce a projection operator e,, whose action 
is to eliminate those terms that contain the factor m in 
the operand. For example, eiUj = (1 — 5ij)uj, etUjUk = 
(1 — 5ij)(l — 5ik)ujUk, etc. As shown below, many-body 
wave functions can be constructed in a systematic way 
by appropriate operation of ej on the yrast-state wave 
function, which is known for < L < N to be Yj, = 

Ei 



<ii<-"<tL<JV 



u h ' ' ' u il [8] ■ The low-lying excitations 
of interacting Bose systems are collective in nature [4]; 
we claim that the main features (i.e., a large energy 



gap and almost linear rotational bands) of the excitation 
spectrum arise from a particular set of coupled collec- 
tive excitations generated by operators P m — Ei^=i u T e i 
(m = 0, 1, • • •). These operators are analogous to Mottel- 
son's multipolar operators Q m cx ^ z™ [4] but differ in 
that Zi is replaced by Ui and in that the projection oper- 
ator ei is incorporated. Both of these modifications are 
essential for identifying an invariant subspace spanned by 
the yrast states and for constructing upon it the desired 
many-body excited states. 

We consider a linear superposition of the states 
M[ m) = P m Y L - m (m = 3, 4, ■■■,£), where M^ m) de- 
scribes an excitation in which one particle carries AM m 
and L — m particles each carry a unit AM, for 



M 



(m) 



E 

l<-<i L - 



L-m+1 

E < 

-1 fc=i 



L-m+l 

1 n * 

1=1 



The remaining TV — L + m — 1 particles carry no AM, and 
thus the condition L < N + m — 1 must be met. 

The transformation law of under the applica- 

tion of V can be obtained as follows. The application of 
27n5(«i — Uj) on m|"^ gives 



(Ui + Uj 



L 

(m) 



(m) , (Ui - Uj'' 2 



Mr + 



eiejM y L _ 2 



Am) 



1m— 1 



-{< + <) 



eiejY L . 



(Uj + Uj) 



m+l 



-(u m Uj+U iU m ) 



eiejY L - m -\. (2) 



Summing both sides of this equation over i and i), 
we obtain for m > 3 

VM { L m) = {(m + 2 3 -" 1 - 4) A + [2 2 - m (m - 1) + 1] 
x (L - m) + 4(1 - 2 1 -™)} M { L m) /2 - (1 - 2 2 ~ m ) 



x(N-L + m)M { L m+1) /2 + 2 2 - m m(L -m + l) 



xM 



(m-l) 



R m+1 + 9(m > 4)R m , 



where V =V — N(N - 1 - L/2) and 



,i-m \ ^ m - 



r(fc) 



Km - 2^ k U m _ k y r rn-kM L _ m+k . 

k=2 y ' 



(3) 



(4) 



Equation (3), together with (4), is the desired recursion 
relation. We see that each mode described by M^ is 

(k) 

coupled with other modes M L with k = 3, 4, • • • , m — 1, 

and k = m + l. Since the coefficient of M^ n+1 ^ in Eq. (3) 
includes a factor N — L + m , the weight of this term 
decreases monotonically with increasing L. By introduc- 
ing the truncation approximation discussed below and 
ignoring contributions from higher-order terms such as 
we obtain a closed set of linear equations that 
can be solved easily even for a large value of m. It turns 
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as higher rotational bands. We also note the coupling 



of (to > 4) to M^' becomes very strong, on the 

order of O(N), near L ~ N. To reproduce the lowest- 

(3) 

lying branch consisting mainly of M L , we have to solve 
the coupled equations for m = 3, 4, • • ■. The solid curve 
in Fig. 3 shows the lowest excitation energy obtained by 
solving the above closed set of equations. It matches the 
G mode rather well. 

In conclusion, we investigated the Nambu-Goldstone 
mode (NGM) of a rotating BEC associated with axisym- 
metry breaking due to vortex nucleation, and identified 
this mode with the lowest-lying envelope comprised of oc- 
tupole branches that are equidistantly displaced by the 
admixture of quadrupole excitations. We found that as 
the angular momentum of the system increases, the NGM 
acquires mass due to its strong coupling to higher rota- 
tional bands. 

M.U. and T.N. acknowledge support by Grants-in- 
Aid for Scientific Research (Grant Nos. 15340129 and 
14740181) by the Ministry of Education, Culture, Sports, 
Science and Technology of Japan. 



FIG. 3: Comparison between the excitation energy of the G 
mode obtained by solving a set of linearized equations for 
m = 3, 4, • • • , 7 (solid curve) and those obtained by the exact 
diagonalization of the Hamiltonian (open circles). The energy 
is measured from the yrast state, which corresponds to the 
horizontal bottom line. 



out that to quantitatively reproduce the G mode labeled 
in Fig. 1, we must take into account to up to a rela- 
tively large value, such as 7, because higher rotational 
bands couple strongly to the G mode, as discussed be- 
low. Writing Eq. (3) explicitly for to = 3,4, and 5, we 
obtain 



vm! 3) = 



--L(L-2)Y L , 



VM 



(4) 



-M 



A3) 



2N + 7L- 24 



M 



(4) 



VM 



(■>) 



|(iV-L + 4)Mf + Il{L-2)Y l , 
■|(L-3)Mf -A( L _2)M< 4 > 



M 



(5) 



107V + 12L- 55 
16 



--(N-L + S)M<*>, 



(5) 



Am) 



We note that the diagonal coefficients of M L m 
Eq. (5), are (2iV + 7i)/8 and (57V + 6£)/8 for m = 4 and 
5, respectively, and that they already reproduce rather 
well the almost linear spectra in Fig. 3. We refer to these 
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